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TECHNICAL PAPER 


THE ESTIMATION ERROR COVARIANCE MATRIX 
FOR THE IDEAL STATE RECONSTRUCTOR 
WITH MEASUREMENT NOISE 


I. INTRODUCTION 


Recently Polites developed a new state reconstructor for deterministic digital control systems 
which he called >he Ideal Sla.e Reconstructor (1SR, 11-51. I. was so named because lor a linear . me 
invariant plant with noise-free measurements, measurable inputs, and known parameters, t y 

reconstructs ^h^sta^te of the plant and does so without requiring any knowledge ot the plant s initial 
nTdd Uon he SR adds no new states or eigenvalues to the system. Nor does it at tect the plant equation 
fo/the systemi n any w ay ; it affects the measurement equation only. However, systems in the real worid 
normally have some measurement noise, nonmeasurable inputs, and parameter uncertain y. 
effect of these things on the ISR needs to be determined. This paper addresses the rst. 

The outline of the paper is as follows. First, a briel review of the ISR is given, m Section I 
Then, an analysis of the ISR with measurement noise is presented, in Section »• N«^an examp ms 
given in Section IV. which illustrates the effect of measurement noise on the ISR. Finally, 
elusions and recommendations for further study are presented, m Section V. 

II A REVIEW OF THE IDEAL STATE RECONSTRUCTOR 


Consider the ubiquitous continuous-time plan, driven by a zero-order-hold with » 
as shown in Figure I matrix G € R nxr ; the standard 

discrete state equations 1 6 ] 


x[(k+l)T] = Ax(kT) + Bu(kT) 


y s (kT) = C s x(kT) 


( 2 ) 


where 


4>(t> =/'' usi-Fr'i e R nxn 


(3) 



( I ) Choose N so that 


N ' n/p 


( 6 ) 


(2) Form the matrix a where 


a 


C, d>(0) 


C,. <*>( - 



£ j^(Np)xn 


C f . <t> I — (N — 1 > — C | 
N 1 


(7) 


and check to be certain it has full rank (i.e., rank n). If so, proceed by letting 


a* = (a'a )' 1 a 1 € R" x(N i” 


( 8 ) 


(3) Choose q so that 


q 5= n - m 


(9) 


(4) Now choose a matrix D_ € R t|xn so that 


Cj 


" C s ‘ 
- D - 


£ \ t.|}\n 


( 10 ) 


has full rank (i.e., rank n). 

(5) Let 


C'r = (C t t C t )'' C t t € R nx<m, M> 
and 


4 


H = D_ a* € R 4X,Np) 


( 12 ) 


(6) Partition H as follows 


H 


LH„ 


H, 


Hn-iI 


to reveal Hj, j — 0,1 N-l. 

(7) Form the matrix p where 


P 


C H [/ <}>(X)dX] G 
0 

-(T/N) 

C, [/ <t><\)d\] G 




)xr 


-<N-I)(T/N> 

C F \f <J>(X)dx] G 


and let 


E_ = Hp € R HXI 


The 1SR is now completely defined for exact state reconstruction. 

A special form of the ISR is obtained when C, is a null matrix (i.e.. m = 0) and q = n - m 
Then, D_ € R nxn and, by virtue of equations (10) and (11), 

C T = D_ € R nxn 


and 


C 


T 


D ' € R nxn 


respectively. 


= n. 

(13) 

(14) 


5 


III. ANALYSIS OF THE IDEAL STATE RECONSTRUCTOR 
WITH MEASUREMENT NOISE 


Now consider the system in Figure 3, which is the system in Figure 2 with 


measurement noise 


S ' S ^ b ; l 'd* n0teS a I"""* corru P ,ed b >' no * se - The measurement noise sources are n s (t) 
and c,.(t) € R I hey are assumed to have the following statistical properties after sampling: 


E ks(kT)] = 0 € R n 


E | c s (kT) (jT)J = R ss 8 ( k - j ) € R 1 


E K( k ^j-)l = 0 € R p , 


E ki (k X> y,. T (j X)] = R|1 8 (k .j) e R rxp 


E k(kT) u,. T (kT - j X), = Rs| 5(j) € R mxp 


E k (kT) .-s'tkT -jl)| = R,, s 8(j) = R sr . T 8(j) e R px 




(15) 


to* dl1 m,e " cr values ol k and j- 11 ls straightforward to show that the output of the ISR in Figure 1 can be 
expressed as 6 * uc 


y r (kT) = y-r'(kT) + Cx 


/ X'X k XL \ 

' H V|.(kT) / 


where y r '(kT) is the output for the noise-free case a 


ind 


Yi(kT) = 


Jd (kT) 


U, ikT - x, 


!’)■ 1 kT - (N-l) Xj 
N 


€ R Np 
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Assume the parameter matrices in the ISR are chosen by the procedure described in Section II and exact 
state reconstruction is achieved in the noise-free case. Then y x '(kT) is given by equation (5). Obviously, 
with measurement noise, y r '(kT) is an estimate of x(kT). Hence, the state estimation error for the ISR can 
be written as 


x 0 (kT) = y T '(kT) - x(kT) = C T * 


| Ts(kT) x 

' h y, (kT) ! 


e R n 


06 ) 


With equation ( 16) and the statistical properties of the measurement noise specified in equations (15), the 
statistical properties of x c (kT) are found to be 


E |x,.(kT)] = 0 


and 


COY [x,(kT), Xe(jT)] 


E [x 0 (kT) x 0 T (jT)] - E [x c (kT)J E [x c T (jT)] 


Rcc S(k-j) , 


for all integer values of k and j, where 


or 


Rss 1 [RsfO ... OJ H t 


Rcc = Cl 


“Rsp 1 " 

1 

1 

Rpp 

0 


H 

0 

! H 

t 





_ 0 

i 

» 

1 

1 

_0 

R,r_ 



Rss 


Rs. H 0 T 


Rcc = or 

Ho R sf t 

N-l 

Z 

(H, R ff H/) 

(C-?) T 



i=0 




(c-fr 


> 



(17) 
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If the performance index J for the ISR is chosen to be the trace of COV[x c (kT), x c (kT)], then 


J = trace (COY |x e (kT), x c (kT)|j = trace (R cc ) 


( 18 ) 


which is invariant with time. 

Consider the special form of the ISR where C s is a null matrix (i.e., m = 0) and q = n-m = n. In 
this case, there are no standard measurements and, hence, no associated measurement noise. Therefore, 
R.ss and R S f are null matrices. With this and equations (12), (14), and (17), 


Rff 0 


Rcc = a* 


(a*) 1 


L 0 Rf, 


FF 


(19) 


By virtue of equations (7), (8), (15), (18), and (19), once the plant is fixed and the measurement noise is 
specified, the performance index J is only a function of N, the number of measurements used in estimat- 
ing the state every T seconds. The sensitivity of J to N will be illustrated by means of an example, pre- 
sented in Section IV. 


IV. AN EXAMPLE 


Consider the system in Figure 4. It is a double integrator plant with measurement noise and the 
special form of the ISR described in Sections II and III. It is assumed that the statistical properties of the 
measurement noise after sampling are given by equations (15) with Rff = L Rss and R S f are null 
matrices. Furthermore, it is assumed that the time interval between state estimates is T = I sec. The 
number of measurements N, used in estimating the state at a given time, is a parameter which will be 
varied to see the affect of N on the estimator performance index J. 

Manipulating the plant in Figure 4 into the format of Figure 3 yields 


F 


0 1 
0 0 


( 20 ) 


G 


0 

1 


9 











and 


C,. = [I 0) , 


where C s is a null matrix. It follows from the definitions of F. G, C, , and C s given in Section II that n 
-■ r = P - •• and m = 0. From equations (3) and (20), 


d>(t) 


1 

0 


( 22 ) 


Proceed now with the design of the ISR using the procedure outlined in Section II. From equation 
(6), the requirement on N is that N 2= 2. From equations (7), (21), and (22), 


I 0 

I -I/N 


a 


-2/N 


I -(N-l)/N 


when T - 1 sec. For N 2= 2, a is seen to have full rank |i.e., rank (a) = 2] as required. Since N is a 
variable, determination ol a* according to equation (8) is best done numerically. From equation (9), the 
requirement on q is that q 2= 2. For the special form of the ISR, q = 2. Choose 


I 0 


D 


0 


so that, by virtue of equation (13), 


Cr 


I 0 
0 I 


which has lull rank as required. 


11 



The remainder of the calculations to define the ISR are straightforward and will be omitted here, 
except for the calculation of the performance index J, defined by equations (18) and ( 19) for the special 
form of the ISR. This was done numerically for various values of N. The results are shown graphically in 
Figure 5. The graph indicates that the more measurements used in estimating the state at a given time, the 
better the estimate, on a statistical basis. It also shows that for N 3* 4, J is essentially proportional to 1/N, 
ala the variance for the average of N independent random variables with the same statistical properties |7J. 



N 

Figure 5. J versus N for the example. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


This paper derived a general expression for the state estimation error covariance matrix for the 
Ideal State Reconstructor (ISR) when the input measurements are corrupted by measurement noise. An 
example was presented which showed that the more measurements used in estimating the state of the 
plant at a given time, the better the estimate, on a statistical basis. Intuitively, one would expect this to be 
true in general. Assuming this is the case, then the ISR is ideally suited for systems where measurements 
become available more frequently than the state needs to be estimated and control law equations need to 
be solved. A good implementation of it would be to have a microprocessor dedicated to multiplying the 

measurements made every T/N seconds by the weighting matrices H r j = 0,1, N-l, and summing 

these results, recursively. After N repetitions, the result could be transferred to a central processor along 
with the measurements made every T seconds. In the central processor, the remainder of the calculations 
to estimate the state of the system could be made and the control law equations solved. Before the control 
law equations are solved, the estimated state from the ISR could be passed through a Kalman filter to 
further improve it. This is an area for future study. Determining the effects of process noise and modeling 
errors on the ISR is also. 
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